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Like concurrence, negativity reflects the monogamous nature of entanglement among multiple 
qubits. Both satisfy a linear inequality involving a party’s entanglement with other parties. We 
show that negativity only saturates the inequality in trivial cases, and we derive a generally tight, 
higher order polynomial inequality. Our methods then shed light on the shareability of negativity 
in larger systems. 


Entanglement has been long recognized as a resource 
for quantum informational tasks [1] and recently, entan¬ 
glement has played a role in manifesting novel phases of 
matter [2]. In each setting, the quantitative theory of 
entanglement aids the analysis. An important result of 
quantitative theory is the strict bounds it places on the 
distributability of entanglement, which can be used to 
deduce various features of quantum systems [3]. 

According to the quantitative theory, establishing cer¬ 
tain physical properties of an entanglement measure 
forces it to identify uniquely with pure state entangle¬ 
ment entropy [4]. In the setting of two pure qubits, the 
entropy is solely determined by an SU(2)xSU(2) poly¬ 
nomial invariant named concurrence. Studying such an 
algebraic invariant alone enabled the solution to the fol¬ 
lowing minimization problem: What is the smallest av¬ 
erage entanglement entropy of an arbitrary mixture? [5] 
(see also [6]). With a mixed state measure of entangle¬ 
ment readily computable, entangled pairs in multiqubit 
systems could be understood, and from there arose the 
notion of entanglement monogamy. 

Entanglement monogamy is the property that the en¬ 
tanglement between Alice’s qubit and Bob’s qubit is 
limited by the entanglement between Alice’s qubit with 
Charlie’s. The property was first captured by concur¬ 
rence in the following monogamy inequality [7], 

^A\BC — ^A\B + ^4|C> (1) 

where C is the concurrence, subscripts label the par¬ 
ties, and the vertical bar indicates the bipartite splitting 
with which to compute the concurrence. In this way, if 
C\\bc ~ ^a\b ~ th en Ca\c ~ hence the personifica¬ 
tion of entanglement. 

To demonstrate the tightness of (1), it is informative to 
develop a pictorial representation. Consider the achiev¬ 
able set of triples, {C\\ G ,C\\ B ,C\\ BC ), computed for all 
pure states of 3 qubits. The equality in (1) then de¬ 
fines a plane in the ambient space, [0, l] 3 . The original 
work [7] found saturation with the W class of states, 

| W) = a|001) + b|010) + c|100), which do, in fact, map 
onto the entire plane in [0, l] 3 . Less well-known is the fact 
that the entire region above the plane can be achieved; 
see Fig. 1. 

The elegance and simplicity of the monogamy inequal- 



FIG. 1: Achievable qubit concurrence 


ity has made it the paragon for entanglement shareabil¬ 
ity. It has since been shown that entanglement of forma¬ 
tion [8], squashed entanglement [9], and entanglement 
negativity [10], all satisfy the same monogamy relation. 

Negativity, in particular, is important for several rea¬ 
sons. It is directly related to PPT states, a peculiar set of 
entangled states that, among other properties [11], can¬ 
not be distilled [12]. Further, negativity provides an al¬ 
ternative measure of mixed state entanglement that has 
the extremely rare property of being computable. In this 
Letter, we show that the ‘linear’ monogamy inequality 
exemplified in (1) fails to describe the achievable set of 
negativity triples and we derive the true boundary. By 
generalizing our approach to qubits, we are then able to 
formulate a conjecture on the strict shareability of nega¬ 
tivity in arbitrarily large quantum systems. 

Negativity is based on the failure of the transpose op¬ 
eration to preserve positivity when acting on subsys¬ 
tems [12]. Transposing a separated system leaves the 
positivity unaffected, so a state with a non-positive par¬ 
tial transpose must be entangled. Negativity is defined as 
twice the sum of the negative eigenvalues of the partially 
transposed state: 

Na\b =N(pab) = 2 ^^ A [, \pab)> ( 2 ) 

n 

where is the density matrix with the A tensor fac¬ 
tor transposed. Although negativity is not obtained by 
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a minimization procedure, it remarkably has many prop¬ 
erties expected of an entanglement measure, e.g., it has 
an interpretation as a distance to a separable state [13]. 
Negativity is sufficient to detect any entanglement in 
2x2, i.e., two qubit systems [14]. Although failing to 
respond faithfully in general, it does bound, e.g., distill¬ 
able entanglement [12] and remains widely used. One can 
straightforwardly verify that for pure 2 x 2 systems, neg¬ 
ativity and concurrence agree, M = C. When extended 
to mixed states, due to the convexity of negativity [ 12 ], 
we have 

C(pab) = inf^p n C(|v4”2)) 

n 

= inf^p„AA(|^i)) 

n 

> inf A r p n | V’is) <V 4"2 1 j 

= N{pab), 

see also [15]. Also, for pure 3-qubit states, Ma\bc = 
Ca\bCi and therefore the same monogamy inequality is 
satisfied, 

•Na\bc'^.Na\b+-Na\C’ (3) 

as noted above. 

In order to find the tighter monogamy inequality, i.e., 
the boundary of the achievable set {J^a\ci^"a\bi^"a\bc)i 
it will be useful to have a parametrization for the 3- 
qubit pure states. Acfn et al. showed how to “rotate out” 
all local unitaries to achieve a canonical form—a tripar¬ 
tite analog to Schmidt decomposition for pure bipartite 
states [16]. One such form is given as 

| S') = d|000) + w|100) + a|101) + 6|110) + c|lll) (4) 

with real a, b, c, d > 0, u> £ C, and the usual normaliza¬ 
tion. The boundary will be computed from “below” by 
maximizing the negativities in this parametrization. 

Qubit negativity necessarily satisfies 0 < J\f < 1. We 
find the one-two party split negativity straight forwardly, 

M"X\ BC = 4(a 2 + b 2 + c 2 )d 2 = 4(1 - d 2 - |w| 2 )d 2 , 

where the last equality follows from the normalization 
constraint. By maximizing with respect to d, we find the 
inequality, 0 < K\bc < (1 - |w| 2 ) 2 . Thus when N\\ B c 
is maximal, the parameter ui will vanish. For N\\ B c not 
maximal, it will be useful to consider for what values of 
ui are and both maximized. 

To calculate \ C , we need the following fact: a par¬ 
tial transpose cannot produce more than (D — l ) 2 neg¬ 
ative eigenvalues for two entangled D-dimensional sys¬ 
tems [17]. Accordingly, for two-qubit states, p 7 ^ has 
no more than one negative eigenvalue. The negativity 


is twice the negative eigenvalue, and thus satisfies the 
quartic polynomial equation, 

0 = det|2p^ +ATa\c^\ 

= —16 a 4 d 4 — 16a 2 c 2 d 4 — 8 a 4 d 2 :r + 8a 2 b 2 d 2 x 

— 8 a 2 c 2 d 2 x — 8 a 2 d 4 x — 8a 2 d 2 \uj\ 2 x + 4 a 2 b 2 x 2 
+ 4 b 2 d 2 x 2 + 4 c 2 d 2 x 2 + 4c 2 |w| 2 :r 2 + 2 a 2 x 3 

+ 2 b 2 x 3 + 2 c 2 x 3 + 2 d 2 x 3 + 2 |w| 2 a : 3 + x A 

— (16d 2 a: + 8x 2 )abc\u\ cos(arg(w)), (5) 

where 1 4 is the 4x4 identity matrix and x = Ma\c- Im¬ 
plicitly differentiating this quartic with respect to arg(w) 
and setting 9x/<9(arg(u;)) = 0 gives 

8 a 6 c|w|a:( 2 d 2 + x) sin(arg(a;)) = 0 , 

so henceforth we restrict u> S R + and drop the absolute 
value; the potential extra minus sign from u> £ K _ does 
not affect the end result. Once again, differentiating the 
quartic in (5) with respect to w and setting dx/du ; = 0 
gives 

x(4 abed 2 + Aa 2 d 2 uj + 2 abex — 2 c 2 wx — ujx 2 ) = 0. ( 6 ) 

The gives one constraint, and along with normalization, 
leaves only three parameters. Since we are after a bound¬ 
ary surface in [ 0 , l] 3 , we will eliminate another variable: 
Using (5) to maximize with respect to c gives: 

8a 2 cd A + 4a 2 cd 2 x + 4 abd 2 u>x — 2cd 2 x 2 

+ 2abu>x 2 — 2 cjjJ 2 x 2 — cx 3 = 0. (7) 

Now we employ a powerful technique from computa¬ 
tional algebraic geometry to perform algebraic elimina¬ 
tion. Finding the minimal generating set, the Grobner 
basis, for the ideal generated by these polynomial con¬ 
straints [18], will give polynomials with the proper vari¬ 
ables eliminated. The Grobner basis for Eqs. 5 (with 
w £ K), 6 and 7 has a single element; setting it to zero 
gives: 

(2 ad — x) ( 2 a 2 + x)(2 ad + x)(2 d 2 + x) 
x ( 2 a 2 d 2 + a 2 x - b 2 x) 2 (Aa 2 d 2 - 2b 2 x - x 2 ) = 0 . ( 8 ) 

Neglecting negative solutions for x leaves three options: 
x = 2 ad, and x = — b 2 + -\/& 4 + 4a 2 d 2 , the former produc¬ 
ing a sub-manifold of the latter with b = 0. The solution 
x = 2a 2 d 2 /(b 2 — a 2 ) also produces a sub-manifold: to 
show this, we now enforce normalization and find one 
more Grobner basis. Eliminating x, u>, c, given the fol¬ 
lowing constraints, x = 2a 2 d 2 /{b 2 — a 2 ), Eqs. 5 (with 
ui £ R), 6 , and normalization, produces, 

a 7 b 3 d 4 (a 2 — b 2 ~ ad){a 2 — b 2 +ad)(—a 2 + b 2 + 2a 2 d 2 ) 2 = 0, 

(9) 
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demonstrating that x = 2a 2 d 2 /(b 2 — a 2 ) = 2 ad is again 
a sub-manifold of solution x = —b 2 + y/b 4 + Aa 2 d 2 . A 
similar analysis on leads to the following triples, 

( N\\ C \ ({b 2 - Vb 4 + 4a 2 d 2 ) 2 \ 

Na\b = (a 2 — \/a 4 + Ab 2 d 2 y ■ (10) 

y^A\BcJ \ 4(a 2 + b 2 )d 2 ) ) 

These triples come from precisely the condition that 
oj = c = 0, leaving states in (4) that are locally equiva¬ 
lent to the W class, the same class that maximize con¬ 
currence. The three components of (10) parametrically 
define the boundary of the achievable set. Together with 
the normalization constraint, we can eliminate the state 
coefficients, turning the parametric surface into an im¬ 
plicit surface. Computing the Grobner basis of the para¬ 
metric polynomials, we again find a single element, which 
set to 0 gives the surface implicitly: 

z 6 - 2z 4 {x 2 - xy + y 2 ) 

+ z 2 [x 4 + y 4 — 2 xy(x(x — 1) + y(y — 1) — 3xy/2 + 2)^ 
+ xy(2y 2 + xy 2 + x 2 y + 2x 2 )(x + y + 2) =0, (11) 

where we identify (x 2 ,y 2 ,z 2 ) = {J^\\ G ^A\b^a\bc)- 
We now show that (11) is the only non-trivial boundary 
of the achievable set. Adding back in the parameter c 
will fill in the rest of the set. The partially transposed 
reduced states are then full rank so that the determinants 
are negative as long as there is entanglement, 

det p T A A c = —a 2 d 4 (a 2 + c 2 ), 

d et PA\B = ~b 2 d 4 (b 2 + c 2 ). 

Since N\\ BC = 4(1 — d 2 )d 2 , fixing z 2 will fix d. 
Thus, on a constant z 2 plane, starting at the bound¬ 
ary (c = 0), which is a curve intersecting the x 2 z 2 — and 
y 2 z 2 —coordinate planes at (z 2 ,0,z 2 ) and (0, z 2 ,z 2 ), re¬ 
spectively, increasing c to \/l — d 2 will smoothly collapse 
the curve into the point (0,0, z 2 ) as a and b , and hence 
the determinants and the negativities, vanish. During 
this collapse, the curve continues to intersect the coor¬ 
dinate planes and traces out the achievable z 2 -plane set. 
As z 2 is arbitrary, all points between the boundary and 
the z 2 -axis can be achieved. See Fig. 2 for the achievable 
negativity set. 

The states whose negativities fill up the region shown 
in Fig. 2 can be understood to have a special form, which 
is useful for generalizing this result to higher dimensional 
tensor factors. From (4) with u = 0, 

|^) = d|000) + a|101) -I- 6|110) + c|lll) 

= a| < I ) ) J 4c|0)B + &| < I ) )ab|0)c + c|GHZ) abc 
+ (d — a — b — c) 1000), 



FIG. 2: Achievable qubit negativity 


where |4>) = |00) + |11) and |GHZ) = |000) + |111). These 
states then generalize to D-dimensional qudits straight¬ 
forwardly via |$) I jj) and I GHZ) I jjj), 


D—l 

I'F) = e?| 000 ) + a \j°j) + b \n°) + c| jjj)- ( 12 ) 

j=i 

The partial transpose of the reduced density operator for 
(12) block diagonalizes to: 


o Ta 

Va\c 


(■D-1XP-2) 

3 =1 




D—l 

© 


j =i 


/0 ad 0 \ 

[ ad b 2 be I 

y 0 be a 2 + c 2 J 


Pa]b ^ as same form with a and b interchanged. The 
similarity with the D = 2 case, particularly the 3x3 
matrix factor, tells us that setting c = 0 maximizes the 
pairwise negativities for this family. The third negativity 
can be computed to be: 


Ua\bc = 2 (D — 1 )d\/ a 2 + b 2 + c 2 

+ (D — 1)(D — 2)(a 2 + b 2 + c 2 ), 

so that again, fixing d will fix Ma\bc- T de same argu¬ 
ment about the determinants applies again, so 

det P Ta c = (-I)LtJ d 2 ((a 2 + c 2 )d 2 ) {D - 1) a 2 ^ D - 1 ^ D - 2 ^ 

detpl A B = (-l)^id 2 {(b 2 + c 2 )d 2 ) iD - 1) b^ D -^ D ~ 2 \ 

on a constant A/)i|BC"Pl ane - ^ nd as before, increasing 
c to \J (1 — d 2 )/(D — 1) will send a and 6 to 0 so the 
determinants, and thus the pairwise negativities, vanish. 
Since this is a natural extension of the achievable region 
we found for D = 2, we conjecture that it is the entire 
achievable set of negativities, {.■^"\\c^a\b^a\bD' ^ or 
D > 2 as well. 
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For the boundary states (c = 0) the parameters can be 
eliminated for the Grobner basis to get the conjectured 
implicit bound; however already in the D = 3 case, the 
polynomial is rather complicated, containing 143 terms. 
It is worth mentioning that naively testing the conjec¬ 
ture numerically is nearly hopeless, since the negativities 
of random states are highly non-uniform throughout the 
achievable set [19]. Nevertheless, testing for perturba¬ 
tions of our boundary has led to no counter-examples. 

An alternative way to derive our boundary states is 
the following. Consider the class of maximally entangled 
states between A and B, with an ancillary qudit, C”, 

D -1 

\^) = (d\00) + Y l b\jj))\0). (13) 

3=1 

These states give the entire line, Ma\b = ASa\bc with 
Ma\c = 0. Furthermore, given the two-qudit swap 
Hamiltonian, 

#SWAP = (J) 
j<k 

where a x k ^ is o x acting in the {j, fc}-subspace (unless j = 
k in which case it is just 1 acting on the {j}-subspace), 
let it act on the B and C qudits in (13): 

e i9H SWAP ^ = de i(, |0 0 0} 

D — l 

+ ^E c °s ( 6 , )|j'jO> + isin(0)| jOj). (14) 
3=1 



FIG. 3: Achieved D = 3 negativity of states in Eq. 14. Red 
indicates the (nonbounding) part of the polynomial surface 
for d < 1/ y/D. 

where the subscripts denote the degeneracy. When d > 
1 /'/D, the remaining, smaller, eigenvalues saturate the 
marginal inequality: 

(D - 1 )(a 2 + b 2 ) <{D- l)a 2 + (D - 1 )b 2 . 

The expressions for the negativities are similar, to leading 
order in D: 


Then the phases can be cleaned up with local unitaries 
to match (12), with c = 0. 

Note that the boundary states (14) produce a surface 
that extends and folds back into the achievable set de¬ 
pending on the value of d , i.e., only maximizing the pair¬ 
wise negativities when d > 1/\]D 1 as seen from para¬ 
metrically plotting the resulting polynomial surface in 
Fig. 3. Recall that negativity for D x D systems has 
bounds 0 < AT < D — 1 [1]. 

The condition on d is related to conditions on the 
marginal eigenvalues: Higuchi found a necessary condi¬ 
tion on the univariate marginal eigenvalues for pure N- 
qudit systems [20]; for three qudits it is, 

E ^ ^ E X n ] + E A i c) - ( 15 ) 

n—1 n—1 n =1 

(P) 

including permutations of the parties, where An ' < 
An + i, n € {1, 1} are the ascending eigenvalues of 

party P’s state. The marginal eigenvalues of (12) with 
c = 0 are 

\^ = {{a 2 + b 2 ) (D _ 1): d 2 } 

A (B) ={(b 2 ) (D _ 1) ,(D-l)a 2 +d 2 } 

X^ 0 ') = {(a 2 )( £) _ 1 ), (D — 1 )b 2 + d 2 }, 


Af A c \ 

9 ( a 

AA A B I 

= D 2 I b 2 

AAA\nc) 

\a 2 + b- 


+ 0(D). (16) 


If our conjecture about achievable negativities is true, 
then in the limit of large dimensions, the monogamy in¬ 
equality simplifies to 


Na\bc > Na\b + Na\Ci 

up to terms of 0(1/D). See Fig. 4 for the achievable 
negativities without squares , (Ma\Ct^'a\Bt^'a\bc)- 
In summary, we’ve seen that although negativity is not 
itself a polynomial invariant in the state coefficients, the 
qubit negativity satisfies a polynomial equation. Simple 
maximization procedures combined with application of 
Grobner basis computations enabled us to derive an ex¬ 
plicit expression for the boundary of the achievable set, 
a polynomial surface. Generalizing the qubit boundary 
states motivated a conjecture for the boundary of the 
achievable set for arbitrary dimensional qudits, a conjec¬ 
ture supported by numerical experimentation. We expect 
our approach to qubits to be relevant in proving the con¬ 
jecture for qutrits, provided the proper 3-qutrit canonical 
form is found. For arbitrary dimensions we suspect an 
intimate connection with marginal eigenvalue constraints 
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FIG. 4: Achievable D —> oo negativity of states (12). 


since our boundary states saturate the eigenvalue bound¬ 
aries. 
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